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Abstract 

We present the program optics that simulates the light response of an ar- 
bitrarily shaped scintillation particle detector. Predicted light responses of 
pure Csl polygonal detectors, plastic scintillator staves, cylindrical plastic 
target scintillators and a Plexiglas light-distribution plate are illustrated. We 
demonstrate how different bulk and surface optical properties of a scintilla- 
tor lead to specific volume and temporal light collection probability distribu- 
tions. High-statistics optics simulations are calibrated against the detector 
responses measured in a custom-made cosmic muon tomography apparatus. 
The presented code can also be used to track particles intersecting complex 
geometrical objects. 
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PROGRAM SUMMARY 

Title of program: optics 
Catalogue identifier: ADxx 

Program Summary URL: http : //www . cpc . cs . qub- 
. ac . uk/cpc . summaries /ADxx 

Program obtainable from: CPC Program Li- 
brary, Queen's University of Belfast, N. Ireland, 
also at http://pibeta.phys.virginia.edu/pub- 
licJitml/optics or the authors. 

Licensing provisions: none 

Computers: Tested on MicroVAX 3100 and DECsta- 
tion 5000/200. Program should be easily portable to 
any UNIX workstation. 

Operating system under which the program has been 
tested: DEC VMS V5.5 system in a batch mode and 
the DEC OSF/1 VI. 3A UNIX environment in both a 
batch mode and with a graphical user interface run- 
ning in a batch or interactive mode 

Programming languages used: FORTRAN and Tk/Tcl 
toolkit 

Program libraries used: CERNLIB programs in 
packlib library and kernlib FORTRAN callable li- 
braries (optional) 

Memory required to execute with typical data set: 
up to 1 Mb with a GUI interface 

Number of bits in a word: 32 (CMOS CPU) or 64 
(Alpha AXP CPU processor) 

Peripherals used: standard input, standard output, 
hard disk, X-terminal (optional), postscript printer 
for graphics (optional, preferably color printer) 

Number of lines in the distributed programs, includ- 
ing test data files and help Hies: 26162 

Program Structure: Code consists of 65 individual 
files containing the subroutines, data files and com- 
mand files. A user can modify or expand the photon 
transport code as well as database files specifying de- 
fault optical properties of the detector surfaces and 
the bulk media. 

Distribution format: uuencoded and compressed tar 



file or, alternatively, compressed tar file 

Keywords: Computer modeling and simulation; 
Monte Carlo simulation of scintillator response; Scin- 
tillation detectors; Computed tomography 

Nature of physical problem: Simulation of the vol- 
ume and temporal light collection probability distri- 
butions given the geometrical shape plus bulk and 
surface optical properties of a scintillation detector. 

Method of solution: The code recognizes cylindri- 
cal, spherical, and parabolical as well as arbitrary 
polygonal scintillator shapes (and optional wrapping 
reflectors) that could couple via light guides or wave- 
length shifters to photosensitive surfaces. The light- 
generating volume can be subdivided into the elemen- 
tary cells. The photons generated within each cell 
arc tracked through the scintillating volume taking 
into account specular, diffuse and rough surface re- 
flections from lateral detector surfaces and wrapping 
reflectors, and the bulk attenuation and scattering 
effects from detector defects 0J^] . 

Restriction on the complexity of the problem: The 
statistical uncertainties of the simulated light collec- 
tion probability distribution are limited by the prac- 
tically tolerable running time (see below). 

Typical running time: The running time depends on 
the number of elementary volume cells chosen and 
the number of scintillating photons generated per 
cell and is therefore problem-dependent. For exam- 
ple, assuming a small-step volume subdivision into 
a 15x15x30 matrix with 6750 elements and aim- 
ing for better than 2% average uncertainty in the 
three-dimensional light nonuniformity function typi- 
cally requires 10 7 photon statistics per cell and run- 
ning time of ~ 24 hours on a 200 MHz computer. 

References: [1] B. K. Wright, Program optics (Uni- 
versity of Virginia, Charlottesville, 1992). 

[2] B. K. Wright, Program tkoptics (University of 
Virginia, Charlottesville, 1994). 
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LONG WRITE-UP 



I. INTRODUCTION 

Scintillation detectors are today widely used in nuclear and particle physics experiments 
for detection of ionizing charged particles as well as photons and neutrons. These detectors 
are also often used in medical instrumentation, process control devices, waste management, 
personal protection and nuclear safeguarding ||. 

In the simplest design, the scintillation detector consists of a scintillator assembly (a solid 
crystal, liquid or gaseous scintillator volume) and an optional lightguide or waveshifting guide 
viewed by a light-sensitive amplifying device such photomultiplier tube (PMT) @. The 
fraction of energy deposited in the scintillator material by penetrating ionizing particles 
is absorbed and re-emitted within several nanoseconds in the form of visible or ultraviolet 
light, so-called fluorescent radiation. The scintillator material itself can be an organic crystal, 
plastics or liquid, inorganic crystal, scintillating gas or glass silicate f|||. Scintillation light 
propagates through the scintillator medium and bounces off its surfaces. A fraction of 
the light eventually reaches the photosensitive area of the light-sensitive device where it is 
converted into an electrical signal. The generated current pulse is typically discriminated 
in a user-designed electronics circuit and integrated for a large number of events to produce 
the pulse-height energy spectrum. 

The main considerations in particle detector design are the total light output and the 
detector timing response, with both factors influencing the energy resolution of the detector. 
The amount of scintillation light reaching the photosensitive surface should be maximized 
and should ideally be independent of the position and linearly dependent on the magni- 
tude of the energy deposition. The timing response is determined by the decay constants 
of scintillating excitation and the size and geometry of the detector itself, as well as the 
characteristics of the photo-sensitive device. The light output of a scintillator primarily 
depends on the conversion efficiency of the deposited ionizing energy to scintillation pho- 
tons. The light collected by the photosensitive surface is also dependent on the surface area, 
the efficiency of the light transport and the scintillator transparency to its own scintillation 
light. The match between spectral response of the photomultiplier, the scintillation light and 
the quantum efficiency of the photocathode are additional controlling factors determining 
photoelectron statistics. 

The magnitude and linearity of the pulse-heights and timing resolution of the scintillator 
signal can be optimized by tuning the light transfer to the photosensitive surface. The 
optimization is executed by the appropriate optical treatment of the detector surfaces and 
by specific choices of shapes and dimensions of the scintillator /light guide assembly and 
different wrapping reflectors. 

The analytical solution of the light transfer and collection process can be obtained easily 
for only a few relatively simple geometries. If the shapes of detectors and lightguides are 
complex or irregular, or if the number of the detector modules is large, the Monte Carlo 
solution of the problem is the most practical and sometimes the only feasible calculational 
technique J7|. 

In order to facilitate the detector design several simple computer programs were devel- 
oped in the past twenty-five years and were used to study photon transport through the 
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elements of the scintillation detectors [fH16|1- These studies focused on the possibilities of 
improving the light collection and the uniformity of light response along the main detector 
axes. Refs. |IB,I7| compared the measurements of pulse-height spectra of the cosmic muons 
and radioactive sources that were moved along the detector surfaces with the simulation re- 
sults. The measurements of the scintillator spatial light output nonuniformities are reported 
in Refs. [p~8| pOf| without comparisons to any matching simulations. 

To our knowledge all published photon transport programs to date have treated the 
properties of rotationally symmetric or regular parallelepipedic detector shapes. The only 
publicly available Monte Carlo photon transport program that is flexible enough to treat a 
larger variety of detector shapes is the CERN Program Library GUIDE7 code ||. 

We initially developed the optics library of subroutines to study the light transport 
in simple cylindrical and trapezoidal scintillator detectors ||21|| . As our need to understand 
the realistic photon propagation in the more complicated detector geometries increased, we 
expanded the code to handle a wider range of detector shapes and included more realistic 
photon-detector boundary interactions. 

In Section we show how the optics subroutines define the detector geometry. We 
explain our method of specifying the physical regions, region boundaries, intersections of 
photon trajectories with the region interfaces and calculation of the normal vectors at the 
points of photon-surface interactions. Section [TIT] touches upon the problem of calculating 
the intersection coordinates between the directed lines and the planar, cylindrical, spherical, 
conical and parabolical surfaces. Section [IV] presents our algorithm for the photon trans- 
port through the detector volume. The overall program structure, including the required 
and optional input parameters and a user-controlled output format is described in Sec. |V]. 
The four progressively more complicated examples of the program applications are given in 
Sec. |VT|. An example of using the optics program in conjunction with the standard GEANT 
code that simulates energy depositions caused by ionizing radiation in particle detectors is 
demonstrated in Sec. |V11| . As an illustration of the physics calculation, we show how to 
simulate the pulse-height spectra of monoenergetic electrons and photons in a pure cesium 
iodide (Csl) calorimeter |22|]. Finally, the instructions for users outlining the required steps 
in the program installation are given in Sec. |V111. 



II. DETECTOR GEOMETRIES 

A. Region Specification 

We define a three-dimensional region by specifying the boundaries of the region and a 
point X0(3) inside the region. A convenient way of storing the region definitions is a FORTRAN 
record structure. While originally a record structure was a VAX FORTRAN extension to the 
ANSI Standard, nowadays almost all FORTRAN compilers provide these extensions The 
record structure of a three-dimensional region is defined in the file REGION_STRUCTURE.TXT: 

C Module REGION.STRUCTURE 



C Record structure for 3-dimensional regions. 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



PARAMETER 



(MAX.SURF = 100) 



STRUCTURE /MATERIAL.STRUCTURE/ 
REAL*4 REFRACTIVE.INDEX 
REAL* 4 ATTENUATION_LENGTH 
REAL* 4 SCATTERING.LENGTH 

END STRUCTURE 

STRUCTURE /REGION.STRUCTURE/ 
REAL*4 X0(3) 



INTEGER* 4 
INTEGER* 4 
INTEGER* 4 
INTEGER* 4 



NPLANE 
NCONE 
NSPHERE 
NPARABOLOID 



RECORD /PLANE.STRUCTURE/ PLANE (MAX.SURF) 

RECORD /CONE.STRUCTURE/ CONE (MAX.SURF) 

RECORD /SPHERE.STRUCTURE/ SPHERE (MAX.SURF) 

RECORD /PARABOLOID.STRUCTURE/ PARABOLOID (MAX.SURF) 



RECORD /MATERIAL.STRUCTURE/ 



MATERIAL 



END STRUCTURE 



I I I I I I I I I 



C End module REGION.STRUCTURE . 

The maximum number of surfaces of a specific type is declared in the user parameter 
MAX_SURF, while the actual number of surfaces of each type bounding a region is speci- 
fied by the integer variables NPLANE, NCONE, NSPHERE or NPARABOLOID. The medium fill- 
ing the region is indexed by an integer variable MATERIAL. For example, the assignments 
REFRACTIVE_INDEX(1)=1.54, REFRACTIVE_INDEX(2)=1 . 00 describe a scintillator medium 
in region 1 and air in region 2, respectively. A database of 17 different materials commonly 
encountered in scintillation detector design is kept in the file MATERIAL_PROPERTIES.TXT. 
Presently defined media properties are the light attenuation length, scattering length, re- 
fractive index, atomic number, atomic weight, and the material volume density. A user can 
add additional materials and media properties or change the default values of the database 
entries. 

Once the record structure of a region is defined, any number of variables having this struc- 
ture format can be declared. Arrays of regions can be defined with the statement RECORD 
/REGION_STRUCTURE/ REGION (100) that creates an array of 100 records, each with the struc- 
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ture format defined by the module REGION_STRUCTURE above. The items in the structure of 
an individual region are referred to as "record, item". The material of the region number 
12 would therefore be referred to as "REGION (12) .MATERIAL" while the y-coordinate of a 
point internal to the region 3 could be accessed through a variable REGION (3) .X0(2). 

Specific two-dimensional and three-dimensional shapes are defined as record structures 
in the file GEOMETRY_STRUCTURES.TXT. This module defines the structures for a point, line, 
circle, plane, cylinder, cone, paraboloid, and sphere, as well as a general polygon structure 
specified by a set of vertices. The comments in the file explain the parameters defining each 
shape. A cylindrical volume is, for example, defined by the following piece of code: 

STRUCTURE /CYLINDER.STRUCTURE/ 

REAL*4 XO (3) ! Point on axis of cylinder. 

REAL*4 A (3) !Unit vector in direction of axis. 

REAL*4 RADIUS ! Radius of cylinder. 

INTEGER* 4 TYPE !Type of surface. 

RECORD /SURFACE_PROPERTIES_STRUCTURE/ PHYSICAL 
END STRUCTURE 

Again, the user can add new geometrical shapes by defining additional structures con- 
forming to the above format. 

The program begins with a call to the subroutine DEFINE_REGIONS where the parameters 
of the problem are either defined or read in from a user data file. An example input format 
for a cylindrical scintillator could be: 

C Read in the geometry parameters (Diameter of the front and back face 
READ (6,*) FR0NT_D1 , BACK.Dl ! of region 1 
READ (6,*) Zl !z-extent of region 1 

C Define material types of regions 

REGION (1) .MATERIAL=2 ! Scintillator 
REGION (2) .MATERIALS !Air 

C Define interior points of regions 
REGION (1) .X0(3)=Zl/2.0 
REGI0N(2) .X0(3)=Zl+2.0 

C Define boundary surfaces of regions 
REGION (1) . NPLANE=2 
REGION (1) .NC0NE=1 

C Detector surfaces 

REGION (1) . PLANE (1) .TYPE=1 
REGION (1) . PLANE (1) .N(3)=-1.0 
REGION (1) . PLANE (1) .X0(3)=0.0 

REGION (1) .PLANE (2) .TYPE=1 
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REGION (1) .PLANE (2) .N(3)=1.0 
REGION (1) .PLANE (2) .X0(3)=Z1 



REGION (1) .CONE(l) .TYPE=1 
REGION (1) .CONE(l) .A(3)=l 

REGION (1) .CONE(l) . SLOPE=(BACK_Dl-FRONT_Dl) / (2 . 0*Z1) 



C Last volume defined contains all the others 
REGION (2) .NSPHERE.=1 
REGION (2) .SPHERE (1) .TYPE=-2 
REGION (2) .SPHERE (1) . RADIUS=2 . 0*Z1 

This piece of code defines a cylindrical scintillator with a zl cm long z-axis and front and 
back face diameters FR0NT_D1 and FR0NT_D2, respectively. The detector is contained inside 
a spherical volume filled with air. The different surface types are described in Sec. pTC. 



B. Polygonal Geometries 

The irregular polygonal volume can be described by a set of vertices stored in the 
POINT_STRUCTURE record: 

STRUCTURE /POLYGON.STRUCTURE/ 

RECORD /POINT.STRUCTURE/ VERTEX(IOO) ! Set of vertices. 
END STRUCTURE 

This format gives the number of vertices of the front and back detector face in the first 
input line, and the x, y, z coordinates of the front (z = 0) and back detector face vertices 
in the second and third lines, respectively. For example, the truncated pentagonal detector 
geometry can be defined with a three-line entry: 

5 

2.45 3.37 0.20 3.55 0.0 0.52 7.11 0.0 0.52 8.21 3.37 0.20 5.33 5.46 0.0 
0.00 6.24 22.02 2.04 0.0 22.61 8.63 0.0 22.61 10.66 6.24 22.02 5.33 10.09 21.66 

The above input file describes a 22 cm long truncated pentagonal pyramid with the front 
(back) face side 3.56 (6.59) cm long. 



C. Surface Specification 

The recognized types of detector surfaces are distinguished by the integer variable TYPE. 
This parameter allows a user to specify the physical properties of the surface. The predefined 
surface types are perfectly transparent, perfectly reflecting, and perfectly absorbing inter- 
faces as well as the "realistic" dielectric surface. These idealized surfaces are distinguished 
in the transport code by the parameter TYPE that is assigned 0, -1, -2 or positive inte- 
ger values. A user should also set the optical properties of a surface through the structure 
SURF ACE_PR0PERT I ES_STRUCTURE: 
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STRUCTURE /SURFACE_PROPERTIES_STRUCTURE/ 
REAL*4 ROUGHNESS ! Number describing roughness. 

REAL*4 REFLECTIVITY ! Fraction of incident light reflected. 

REAL*4 DIFFUSE_FRACTION ! Diffuse.ref lection/total.ref lection . 

END STRUCTURE 

The REFLECTIVITY values could range from 0.0 (a perfect absorber) to 1.0 (a perfect reflec- 
tor). The DIFFUSE_FRACTION coefficient determines the ratio between the specular reflection 
(DIFFUSE_FRACTION=0.0) and the diffuse scattering processes (DIFFUSE_FRACTION=1.0, see 
Sec. [IV B|) . The ROUGHNESS parameter describes the deviation of the detector sides from an 



ideal flat surface behavior (Sec. |1V C|) . 

It is important to note that a pair of adjacent regions that are not nested is always 
separated by the two surfaces. This arrangement allows for a boundary between two regions 
to be of a different type, depending on the region side from which a photon is approaching 
the interface. An example is a "two-way mirror" between two regions that could be created 
by making the surface of one region perfectly reflecting and the adjacent surface of the 
neighboring region completely transparent. 



The four different detector geometries are described as examples in more detail in Sec. [V] 
and shown in Fig. |IJ. The parameters of four geometries are specified in Table [TJ These 
shapes represent the real detectors whose light responses we studied both in Monte Carlo 
simulations and in calibration measurements using cosmic muons, monoenergetic tagged 
positron and photon beams, and laser pulses. The geometries are as varied as irregular 
truncated pyramids made of pure cesium iodide, plastic veto scintillator staves, a cylindrical 
plastic target, and a rectangular Plexiglas monitoring plate, respectively. 



III. INTERSECTION OF DIRECTED LINES WITH OPTICS SURFACES 

A straight line in three-dimensional space can be represented by the parametric equation: 

X = X 0L + rV, (3.1) 

where X l is a point on the line, V is a unit vector in the direction of the line, and the 
parameter r specifies the distance from a fixed point Xql to a point X. 

The parametric equations of the plane, spherical, and cylindrical surfaces are given by 

n • (X - X 0P ) = 0, plane, (3.2) 
\X - X os \ = R, sphere, (3.3) 
\X - X oc — [A • (X — X oc )]\ = 0, cylinder, (3.4) 

where n is a normal vector to the plane, X s and R are the center and the radius of the 
sphere, and X c, A are the point on the cylinder axis and a unit vector in the direction of 
the detector z axis, respectively. 

Intersections of a straight line with these surfaces are found in, for exam- 
ple, Ref. |!J. The formulas for the intersection coordinates are implemented in 
the optics subroutine SUBROUTINE FIND _LINE_REGI0N_INTERSECTI0N (LINE , I_REGI0N , 
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XI , X2 , INTERSECT , REGION , N_REGIONS) . The relevant regions should first be defined with 
calls to SUBROUTINE DEFINE_REGIONS_CRYSTAL (REGION , N_REGI0NS) . 

The LINE variable is specified by the three-dimensional array LINE.XO and by the three 
cosines of directions LINE.V that could be calculated from the polar and azimuthal angles 

LINE. V(l) = sin 6 cos <p, (3.5) 
LINE.V(2) = sin 9 sin <fr, (3.6) 
LINE.V(3) = cos#. (3.7) 

The variable I_REGI0N is an index of a region that takes values from 1 to N_REGI0NS, while 
the REGION variable is the FORTRAN record structure described in Sec. Ill A. 



INTERSECT is a logical variable that is .TRUE, if at least one intersection point is found 
and .FALSE, if there are no intersections between the line and user defined surfaces. The 
arguments XI and X2 are three element arrays that contain the calculated coordinates of the 
intersection points. 

In Fig. ^ we show the intersections of reconstructed cosmic muon trajectories with the sur- 
faces of a cesium iodide scintillation crystal. The crystal's geometrical shape is a truncated 
irregular hexagon labeled HEX- A (Fig. |].i). Fig. ^illustrates the cosmic muon trajectories 
intersecting three vertically stacked rectangular scintillator bars (Fig. |l|.ii). In both cases the 
intersection coordinates have been calculated using the FIND_LINE_REGION_INTERSECTION 
subroutine. 



IV. PHOTON TRACKING 

A. Dielectric Reflection 

The reflection and refraction of light from a perfect dielectric interface can described 
by the laws of geometrical optics if the scintillation photons are unpolarized. If the angles 
between the directions of the incident, reflected and transmitted photon with respect to the 
surface normal are denoted with 9^ 8 r and 9 t , respectively, the law of reflection states that 

r = vr - 9i, (4.1) 

while the law of refraction (Snell's law) on the other hand requires that 

n\ sin Oi = n 2 sin 9 t , (4-2) 

where n\ and n 2 are the indices of refraction of the media in regions 1 and 2, respectively. 

The user can change the program so that the full polarization of the propagating photon 
is taken into account. The reflection and transmission phenomena would then be described 
by the Fresnel formulas (see, for example, |§|24|)- 
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B. Diffuse Reflection 



A perfectly diffuse reflector is defined by the uniform scattering probability P(9) into a 
unit solid angle dQ, where 9 denotes the angle between the direction of a reflected photon 
and the normal vector to the surface. The angular distribution of the scattered radiation is 
then isotropic. If the reflecting surface is a plane, the resulting distribution is often called 
the Lambert's (or cosine) law [25|: 

I(9) = I cos9. (4.3) 



C. Rough Surfaces 

A rough surface can be modeled by introducing a set of small plane surfaces tangent to 
the original surface at suitable points. This description goes under the name of the facet 
model and is described in more detail in Ref. Jl6| where the model predictions are also 
compared with the experimental measurements. 

In the simplest version of the facet model a rough surface is characterized by a single 
roughness parameter R. This parameter is zero for smooth scintillator surfaces with a perfect 
finish. For real detectors with uneven boundaries R value is a positive number between 0.0 
and 1.0. Polished rough surfaces with height variations of the order of a micron have the 
roughness coefficient on the order of 0.1. 

The subroutine PERTURB _UNIT_VECT0R (V_IN , ROUGHNESS , V_0UT) calculates a random 
perturbation to the direction of an incident photon unit vector V_IN following the reflection 
from a rough surface. The scalar product of the perturbed vector V_0UT and the unperturbed 
photon direction cos 9 is equal to 1/ ^{1 + 5 2 ), where 5 has a Gaussian distribution, centered 
at zero, with the width given by the ROUGHNESS parameter R. 



D. Absorption and Rescattering 

The media properties are described by two parameters, SCATTERING_LENGTH and 
ATTENUATION_LENGTH. The code simplifies photon transport by assuming that the scatter- 
ing length and the absorption length do not depend on the wavelength of the propagating 
light. The dispersion effects are neglected. If a user wants to consider dispersion effects, the 
changes should be made in the file optics.f. 



V. PROGRAM STRUCTURE 
A. Program Input 

The program starts with the subroutine DEFINE_REGIONS that reads a data file specifying 
the parameters of the detector. The call to that subroutine defines the record structures 
in the array REGION. A starting position of a scintillation photon is accepted next. The 
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coordinate values can be read in from the data file PHOTONS .DAT or, alternatively, a user can 
select one of the predefined distributions: 

1. a single point source, fixed x,y,z; 

2. a through-going particle, range of x, fixed z, y = 0; 

3. a range of light sources on the detector axis, range of z, x = y = 0; 

4. a uniform light distribution throughout detector volume, range of x, y, z\ 

5. a uniformly distributed light sources in the center plane, range of x,z, y — 0; 

6. an aimed "pencil" of scintillating radiation, fixed x,y,z, and fixed directional cosines 

^xi Uyi "zi 

7. a uniform distribution in the first half (first half-sector) of the detector, range of x, y, z. 

All input parameters defining the detector and the scintillating light distribution, as well 
as the desired output format can also be entered from the tkoptics window if one is using 
the GUI version of the program (see Fig. |]). 

B. Physics Subroutines 

The initial region number is found next by calling the routine FIND_REGION_NUMBER. The 
program then steps through the following levels: 

1. Find the next intersection of a photon trajectory with the predefined surface. 

2. Test to see if a photon is absorbed or rescattered in transit. 

2.1 If absorption or rescattering occurred, look at the current region number. 

2.1.1 If a photon is in the wavelength shifter, find the point at which it is absorbed, 
generate a new random direction, and continue. 

2.1.2 If a photon is not in the wavelength shifter, count it as absorbed and get a new 
particle. 

2.2 If a photon is not absorbed, continue. 

3. Propagate a photon to the next intersection. 

4. Look into the next region to find its region number. 

4.1 If the next region number is the same as the current region number, continue the 
photon transport and return to (1). 

4.2 If we are dealing with a different region, continue to (5). 

5. Look at the type of the surface a photon encounters. 

5.1 If the surface is transparent, continue to the next intersection. 

5.2 If the surface is reflecting, stay in the current region while changing the photon's 
direction. 
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5.3 If the surface is absorbing, stop the photon propagation, count it as absorbed, and 
generate a new photon ab initio and return to step (1). 

5.4 When an interface is a real surface, use the Snell's law to determine if the photon 
undergoes reflection or refraction, change the direction appropriately, and step into the next 
region if a photon is refracted. 

5.5 For the watch points, call the WATCHPOINT subroutine, then do (5.1). 

5.6 For the detectors, call the DETECTOR subroutine, count a photon as detected, then do 
(5.3). 

6. Return to (1). 

This algorithm is represented by a box diagram in Fig. |5|. 

The following is the list of subroutines used in the photon propagation: 

SUBROUTINE INITIALIZE ( REGION, N.REGIONS ) 
SUBROUTINE WRITE.STATS ( REGION, N.REGIONS ) 
SUBROUTINE GENERATE.PHOTON ( LINE, REGION, N.REGIONS ) 
SUBROUTINE L0G_PH0T0N_GENERATI0N ( LINE ) 

SUBROUTINE LOG.INTERSECTION ( LINE, REGION, N.REGIONS, N_CURRENT, 

R, X, OUTWARD.NORMAL, CLASS, ISURF, TYPE ) 
SUBROUTINE L0G_REFLECTI0N ( LINE, REGION, N.REGIONS, N_ CURRENT , 

R, X, OUTWARD.NORMAL, CLASS, ISURF, TYPE ) 
SUBROUTINE L0G_REFRACTI0N ( LINE, REGION, N.REGIONS, N_ CURRENT , 

R, X, OUTWARD.NORMAL, CLASS, ISURF, TYPE ) 
SUBROUTINE L0G_ABS0RPTI0N ( LINE, REGION, N.REGIONS, N_ CURRENT , 

R, X, OUTWARD.NORMAL, CLASS, ISURF, TYPE ) 
SUBROUTINE LOG.ATTENUATION ( REGION, N.CURRENT, LINE, X ) 
SUBROUTINE DEFINE_REGION_HALF_SECTOR ( FILENAME, 
SUBROUTINE DEFINE.REGIONS ( REGION, N.REGIONS ) 
SUBROUTINE DEFINE_REGIONS_PYRAMID ( REGION, N.REGIONS ) 
SUBROUTINE DEFINE_REGIONS_CONE ( REGION, N.REGIONS ) 
SUBROUTINE DETECTOR ( TYPE, LINE, RANGE, TIME ) 
SUBROUTINE WATCHPOINT ( TYPE, LINE, RANGE, TIME ) 
SUBROUTINE TERMINATE ( REGION, N.REGIONS ) 

The functions these subroutines perform are described in Sec. |IV[ 



C. Program Input 

The listing of the program input for the graphical user interface (GUI) version of the 
code is printed by default at the program termination: 

Enter geometry type: 1 

= Truncated cone 

1 = Pi-Beta Csl crystal 
Enter pyramid type: 

1 = Pentagon 
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2 = Hexagon-A 

3 = Hexagon-B 

4 = Hexagon-C 

5 = Hexagon-D 

6 = Half Hexagon-D (#1) 

7 = Half Hexagon-D (#2) 

= User defined (will prompt for file name) 
Enter file name: /f s24/users/emil/tomo/tkoptics/hhl-d-mc .dat 
Enter refractive index of pyramid material: 2.10 
Enter attenuation length of pyramid material (cm): 200.0 
Enter scattering length of pyramid material (cm): 200.0 
Enter surface type for sides of pyramid: -1 
1 = normal dielectric 

-1 = imperfect specular reflector 

-2 = perfect absorber 

-3 = imperfect diffuse reflector 
Enter reflectivity of this surface (0-1): 0.975 
Enter roughness of this surface (0-1): 0.0 
Enter diffuse fraction for this surface (0-1): 0.0 
Enter surface type for front face of pyramid: -1 
1 = normal dielectric 

-1 = imperfect specular reflector 

-2 = perfect absorber 

-3 = imperfect diffuse reflector 
Enter reflectivity of this surface (0-1): 0.975 
Enter roughness of this surface (0-1): 0.0 
Enter diffuse fraction for this surface (0-1): 0.0 
Enter refractive index of air gap material: 1.000273 
Enter attenuation length of air gap material (cm): 1.0E+3 
Enter scattering length of air gap material (cm): 1.0E+3 
Enter surface type for pyramid wrapper: -3 
1 = normal dielectric 

-1 = imperfect specular reflector 

-2 = perfect absorber 

-3 = imperfect diffuse reflector 
Enter reflectivity of this surface (0-1): 0.9 
Enter roughness of this surface (0-1): 0.0 
Enter diffuse fraction for this surface (0-1): 1.0 
Enter surface type for front face of wrapper: -3 
1 = normal dielectric 

-1 = imperfect specular reflector 

-2 = perfect absorber 

-3 = imperfect diffuse reflector 
Enter reflectivity of this surface (0-1): 0.2 
Enter roughness of this surface (0-1): 0.0 
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Enter diffuse fraction for this surface (0-1): 1.0 
Enter refractive index of optical joint material: 1.58 
Enter attenuation length of optical joint material (cm): 200.0 
Enter scattering length of optical joint material (cm): 200.0 
Enter refractive index of PMT window material: 1.458 
Enter attenuation length of PMT window material (cm): 100.0 
Enter scattering length of PMT window material (cm): 100.0 
********* Xfront = 6.790000 2.865000 22.19750 

********* AXIS = 0.0000000E+00 0.1182974 0.9929782 
********* VECTOR = 8.4076812E-03 -5 . 3805402E-03 -6 . 0056243E-04 
********* VECTOR = 0.0000000E+00 1 . 1829736E-02 9 . 9297822E-02 
********* xfront = 6.790000 2.876830 22.29680 

Enter output data type: 3 

= No output 

1 = Endpoint data 

2 = Tracks 

3 = Statistics 

Enter output file name: /f s24/users/emil/tomo/tkoptics/tkoptics . rz 

Enter numbers of bins: (NXBINS NYBINS NZBINS) 14 12 24 

Enter X limits of container volume: (MIN MAX) 14 

Enter Y limits of container volume: (MIN MAX) 12 

Enter Z limits of container volume: (MIN MAX) 24 

Enter starting distribution: 3 

= Single point (Fixed X,Y,Z) 

1 = Through-going particle (Range of X; fixed Z; Y=0) 

2 = Range of points along axis (Range of Z; X=Y=0) 

3 = Uniform throughout volume (Range of X,Y,Z) 

4 = Uniform in center plane (Range of X,Z; Y=0) 

5 = Aimed (Fixed X,Y,Z; Fixed VX,VY,VZ) 

6 = Uniform in first half-sector (Range of X,Y,Z) 
Enter number of photons: 5000000 

Number of photons processed = 5000000 

Number of photons attenuated = 1467024 
Number of photons absorbed = 2197182 
Number of photons overbounced = 169244 

Number of photons detected = 1166550 
( 489568 started backward ) 

( 676982 started forward ) 



This output is an example of the optics simulation of a half-hexagonal Csl detector, 
Fig. [TJi. 
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D. Program Output 



The user can choose three different output formats. In the ENDPOINTS format for each 
detected scintillation photon a record is written to an ASCII output file documenting: 

1. xq, yo, z , the starting coordinates of a generated photon; 

2- v xo , v yo , v zo , the components of the initial directional vector of the photon; 

3. X\,yx,zi, a point at which the photon hit the photosensitive surface; 

4. v xl ,v yi ,v Zl , the components of the photon's direction vector when it hit the photosen- 
sitive surface; 

5. t, the time elapsed from the photon creation to the moment of its detection; 

6. n, the number of times a photon was reflected from the detector surfaces. 

The program terminates by printing the photon transport statistics summary. The sum- 
mary includes the number of detected photons, as well as fractions of attenuated, absorbed 
and "overbounced" photons (see the tail of Sec. [VQ listing). The limit to the number of 



photon bounces is set by a user. The default limit is 200 reflections. Events exceeding that 
limit are counted as overbounced and the propagation resumes with a new photon at the 
starting position. 

The CERNLIB PAW macro program titled tkpawread.kumac can be used to read in the 
ASCII file produced by the optics code and make an HB00K4 RZ file containing a sum- 
mary Ntuple ||26|| . The usage syntax is exec tkpawread infile outfile. The structure of the 



generated Ntuple is given in an example below: 

* NTUPLE ID= 10 ENTRIES= 233209 Endpoint Data * 



* 


Var numb 


* 


Name 


* 




Lower 


* 




Upper 


* 




* 


1 


* 


xstart 


* 





180037E+00 


* 





138114E+02 


* 


* 


2 


* 


ystart 


* 





146866E-03 


* 





585641E+01 


* 


* 


3 


* 


z st art 


* 





561714E-02 


* 





225345E+02 


* 


* 


4 


* 


vxstart 


* 




999965E+00 


* 





999996E+00 


* 


* 


5 


* 


vystart 


* 




999988E+00 


* 





999991E+00 


* 


* 


6 


* 


vzstart 


* 




999990E+00 


* 





999998E+00 


* 


* 


7 


* 


xstop 


* 





314337E+01 


* 





115353E+02 


* 


* 


8 


* 


ystop 


* 





238158E+00 


* 





541021E+01 


* 


* 


9 


* 


zstop 


* 





220957E+02 


* 





227119E+02 


* 


* 


10 


* 


vxstop 


* 




999665E+00 


* 





998784E+00 


* 


* 


11 


* 


vystop 


* 




990139E+00 


* 





999993E+00 


* 


* 


12 


* 


vzstop 


* 




994538E+00 


* 





999994E+00 


* 


* 


13 


* 


time 


* 





976813E-02 


* 





302494E+02 


* 


* 


14 


* 


bounces 


* 





000000E+00 


* 





109000E+03 


* 
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The variable names correspond to the ASCII file variables. The number of Ntuple entries 
in this case is equal to the number of generated photons. 

If a user wants to calculate the three-dimensional light collection probability function 
the STATISTICS output format is more appropriate. The variables written to an ASCII file 
and into the resulting Ntuple are: 

1. x,y,z, the coordinates of the centers of the elementary cells into which the detector 
volume is subdivided; 

2. counts, the number of detected photons originating from a given elementary cell; 

3. edge, the flag indicating if an elementary cell is fully contained inside the detector 
volume (edge=l), or is partly outside the detector (edge=0); 

4. thrown, the number of the generated photons having the origins within a given ele- 
mentary cell; 

5. time, the average time-of-flight of the photon (time elapsed between the generation 
and detection moment). 

* NTUPLE ID= 10 ENTRIES= 4032 Statistics * 



* Var numb * Name * Lower * Upper * 



* 


1 


* 


X 


* 


0. 


.500000E+00 


* 


0. 


. 135000E+02 


* 


* 


2 


* 


y 


* 


0. 


500000E+00 


* 


0, 


.115000E+02 


* 


* 


3 


* 


z 


* 


0. 


.500000E+00 


* 


0, 


.235000E+02 


* 


* 


4 


* 


counts 


* 


0. 


000000E+00 


* 


0. 


.957900E+04 


* 


* 


5 


* 


edge 


* 


0. 


•000000E+00 


* 


0. 


, 100000E+01 


* 


* 


6 


* 


thrown 


* 


0. 


•000000E+00 


* 


0. 


.407280E+05 


* 


* 


7 


* 


time 


* 


0. 


000000E+00 


* 


0. 


.463954E+01 


* 



The number of the Ntuple entries represents the number of elementary cells, and in the 
case above corresponds to the 14x12x24 cm 3 subdivision. The simulated detector was 
irregular trapezoidal pyramid shown in Fig. |T].i. The resulting histograms of the number of 
photon reflections as well as the distribution of the photon time-of-flight times are shown in 
Fig. [| Converting the later histogram into a cumulative representation we get the intrinsic 
contribution of the photon arrival times to the timing lineshape of the detector system, 
Fig-0. 

Finally, the information on the individual photon trajectories could be captured with 
the TRACK output format: 

* Ntuple ID = 10 Entries = 53424 Track Data 

* Var numb * Type * Packing * Range * Block * Name * 
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* 


1 


* 


1*4 


* 


* 


* 


2 


* 


1*4 


* 


* 


* 


3 


* 


R*4 


* 


* 


* 


4 


* 


R*4 


* 


* 


* 


5 


* 


R*4 


* 


* 


* 


6 


* 


R*4 


* 


* 


* 


7 


* 


R*4 


* 


* 


* 


8 


* 


R*4 


* 


* 



* 


TRACK 


* 


ntrack 


* 


* 


TRACK 


* 


itrack(ntrack) 


* 


* 


TRACK 


* 


xtrack (ntrack) 


* 


* 


TRACK 


* 


ytrack (ntrack) 


* 


* 


TRACK 


* 


ztrack (ntrack) 


* 


* 


TRACK 


* 


vxtrack (ntrack) 


* 


* 


TRACK 


* 


vytrack (ntrack) 


* 


* 


TRACK 


* 


vztrack (ntrack) 


* 



* Block * Entries * Unpacked * Packed * Packing Factor * 

* TRACK * 53424 * 5604 * Var. * Variable * 

* Total * * 5604 * Var. * Variable * 

* Blocks = 1 Variables = 8 Max. Columns = 1401 * 

The track number, starting and endpoint coordinates and the initial track direction of a 
photon are recorded in this file. Figure [8] illustrates the uniform distribution of the detected 
photons on the surface of a photocathode. The simulated detector was again the HEX-D1/2 
shape shown in Fig. |I].i. 



VI. PROGRAM EXAMPLES 

The coordinate system is always defined with the z axis pointing along the detector 
axis and the origin is placed at the front detector face. If the detector shape is polygonal, 
one vertex of the front detector face should always lie on the positive x axis. All detector 
surfaces in examples below are normal dielectric interfaces. Some surfaces are wrapped in 
the specularly reflecting layer of aluminized Mylar or multiple layers of Teflon sheet. User- 
defined geometrical shapes are included in an encompassing sphere filled with air. The 
rectangular subvolume of the sphere is divided into a number of elementary cells with side 
lengths L, usually ~1 cm. 



A. Irregular Pentagonal, Hexagonal and Trapezoidal Scintillators 

The coordinate system is as follows: the z direction points along the axis of the crystal, 
with z = at the smaller (front) face, and z = L at the larger (back) face. One vertex of 
the front face and one of the back face lie on the positive x axis. The coordinates of the 
vertices for each crystal type are stored in the text data files. The nine different scintillator 
shapes are examined: four irregular hexagonal truncated pyramids (we label them HEX-A, 
HEX-B, HEX-C, and HEX-D), one regular pentagonal (PENT) and two irregular half- 
hexagonal truncated pyramids (HEX-D 1 and HEX-D2), and two trapezohedrons (VET-1, 
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VET-2) |22|. The volumes of studied Csl crystals varied from 797 cm 3 (HEX-D1/2) to 
1718 cm 3 (HEX-C), see Fig. |l[i. The size of the elementary cell was 1cm 3 . 

The geometry of the simulation is simple. The Csl scintillator crystal is wrapped in a 
layer of aluminized Mylar and attached via a 0.2 mm thick optical joint to a 1mm thick 
photomultiplier window. Behind this window is a photocathode. The diameter of the 
photosensitive surface is 46 mm (HEX-D and VET) or 67 mm inches (PENT and HEX) . The 
entire detector is enclosed in a large absorbing sphere filled with air. The program output 
is written to a binary RZ file. The results of the calculation are the three-dimensional light 
collection probability functions illustrated on Figs. ^ and |TD[ The percentage probability 
that photons generated within two-dimensional x-z bins will generate a photoelectron on 
the PMT cathode is shown on top panels. The bottom panels are tomography measurements 
of the corresponding light responses for two representative HEX- A and HEX-D 1 detectors. 



B. Rectangular Plastic Veto 

The cylindrical plastic hodoscope of the PIBETA detector [22] used in charged particle 
tracking consists of 20 rectangular staves, Fig |l].ii. The individual sections are optically 
isolated from each other by one layer of aluminized Mylar wrapping and are viewed by one 
inch phototubes coupled via lightguides on both ends. The simulated light response of the 
individual PMTs as a function of the axial position along the plastic stave is an example of 
the calculation of physical interest. 



C. Cylindrical Plastic Target 



A plastic scintillator cylindrical target is used in the center of the PIBETA detector |22 
to stop a pion beam (Fig. |l|.iii). The positrons from pion and muon decays emanating 
from the target center are detected subsequently in the plastic veto hodoscope and the Csl 
calorimeter. A typical optics calculation of the photon transport predicts the light response 
of the target as a function of the polar angle of the exiting positron. 



D. Light-Distribution Plexiglas Plate 



The lead glass calorimeter monitor for the RadPhi experiment E94-026 at JLAB [27| is 
a laser-based calorimeter monitoring system. The nitrogen laser excites a cylindrical plastic 
scintillator and a generated light is distributed via six fiber optics cables to the 1.2 cm thick 
Plexiglas sheet [^8] facing the lead glass wall in a light-tight enclosure (Fig|l|.iv). The sheet is 



oversized at 152cmxl52cm to minimize the optical nonuniformities. The two-dimensional 
response of the smaller prototype system predicted with the optics code, as viewed through 
lead glass modules is displayed on Fig. O. The maximum predicted variation in the light 
output of the plate as a function of the position of an individual lead glass detector is ~30 %. 
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VII. ENERGY LINESHAPE MONTE CARLO 



A. Interpolation of 3D Light Nonuniformity Function 

The program optics can be used to calculate the light collection probability func- 
tion f(iL,jL,kL), where k are integers and L the step size, on a user-defined three- 
dimensional grid. 

Starting from this three-dimensional probability function the code can calculate for an ar- 
bitrary point source of light (x, y, z) inside the detector volume the light fraction reaching the 
photosensitive device. This is accomplished by means of a series of simple one-dimensional 
interpolations. Our method of choice is the bilinear interpolation on an elementary grid 
square. The accuracy will be improved if one uses not only the light output values but also 
the gradients and the cross derivatives of the light output probability changes. These higher 
order methods that result in a smoother behavior of our function go under the names of the 
bi-cubic interpolations or bi-cubic splines. Numerical recipes for these methods are given as 
FORTRAN subroutines and C functions in Ref. p9|j30| . 



B. GEANT Simulation of Electromagnetic Showers 

The three-dimensional distribution of energy deposited in the detector volumes by ion- 



izing charged particles was calculated by the standard code GEANT [31|. We considered mo- 
noenergetic positrons and photons of 70 MeV total energy incident on the PIBETA calorime- 
ter made of pure Csl pyramids. In the user-written subroutine the energy deposited in every 
step was multiplied by the normalized value of the light collection nonuniformity function to 
yield the detected energy deposition. Sums are kept of both the total deposited energy and 
total detected light. Examples of predicted energy deposition lineshapes in the calorime- 
ter are shown in Fig. [12]. Four different simulations are plotted. The top panel shows a 
pulse-height spectrum for 10 5 simulated photons incident on the Csl calorimeter from its 
center. The responses of both idealized uniform detector and the detector characterized 
with the optics-calculated nonuniformity functions are displayed. The bottom panel shows 
the difference in the response for the 70 MeV positron showers. 



VIII. INSTALLATION AND TESTING 



The computer code described in this paper is distributed in the form of uuencoded 
compressed tar archives named optics. uu (229 kb) and tkoptics.uu (2.4 Mb). These 
two files can be obtained via transfers from the URL http : //pibeta . phys . Virginia . edu- 
/public_html/optics. Same files can also be obtained directly from the authors. The 
UNIX utilities for the tar archive processing can be copied from the official GNU site at 
http : //www. gnu. ai .mit . edu. 

Installation of the optics software on a computer running the UNIX operating system 
is straightforward. Execution of the csh script uudecode taking optics or tkoptics as 
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a parameter creates a compressed tar archive file with an extension tar.gz. The archive 
can be unpacked with the command gunzip. The unpacked optics library files should 
be placed under the /common/local/optics subdirectory. The extracted tkoptics files 
should be copied to the optics subdirectory on the user's disk. If the user prefers the GUI 
version of the program, the installation of CERNLIB libraries, as well as of the Tk and Tel 
libraries |32| is required. The publicly available distributions can be accessed, for example, 
at ftp://ftp.scriptics.com/pub/tcl/. The parameter file .tkoptics defining the input 
parameters echoed in the default optics window should be placed in the login directory. 

The Tk and Tel software should be installed under /common/local/tk and /common- 
/local/tcl areas, respectively. The most recent CERN Program Library release containing 
subroutines used in the tkoptics program is available at the site ftp : // asisf tp . cern . ch. 

The working area can be now named by the project's title and should contain the required 
input files, like the optics.dat and tkoptics.dat. To verify that the program executes 
properly the user should remake the executable code by compiling and linking the source 
code with the make optics command . The Makefile file is included in the distribution. 
The program can then be started with the tkoptics command invoking a GUI interface 
window or it can be run in batch mode with the optics command. 
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TABLES 



TABLE I. Sample ranges of input optical parameters for optics simulations of various de- 
tectors. We simulated the response of the hexagonal and half-hexagonal pure Csl scintillators, a 
plastic scintillator stave, a cylindrical active plastic target, and a Plexiglas light distribution plate. 



Parameter Csl Detector Plastic Hodoscope Plastic Target 1.3 cm Plexiglas 

HEX-A,Dl/2 4.1x0.3x59.8cm 3 04.0x5.0cm 2 100x185cm 2 



Matrix Size 


12x14x24 


4x1x60 


10x10x10 


1x3x10 


Detector Surface Type 


specular 


specular 


specular 


specular 


Detector Surface Roughness 


0.0-1.0 


0.0-1.0 


0.0-1.0 


0.0-0.5 


Wrapper Surface Type 


Teflon 


Al Mylar 


Al Mylar 


Al Mylar 


Photocathode Size (cm) 


04.5/07.0 


01.8 


01.0 


03.0 


Photocath. QE (%)/Timing (m 


20/0.5 


28/0.4 


20/0.5 


13/2 


Bulk Attenuation Length (cm) 


50-300 


100-400 


100-400 


100-400 


Bulk Scattering Length (cm) 


50-300 


100-400 


100-400 


100-400 


Index of Refraction (blue light) 


2.08 


1.58 


1.58 


1.58 


Specular Reflectivity 


0.7-1.0 


0.7-1.0 


0.7-1.0 


0.7-1.0 


Wrapper Reflectivity 


0.5-1.0 


0.5-1.0 


0.5-1.0 


0.5-1.0 
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FIGURES 



FIG. 1. Some of the particle detectors and light distribution systems studied with the optics 
code: (i) HEX-A pure Csl detector, (ii) plastic scintillator veto stave, (hi) cylindrical active plastic 
target, and (iv) Plexiglas light-distribution monitoring plate. 

FIG. 2. The HEX-A scintillator volume intersected by the cosmic muons tracked using a triple 
drift chamber system. The physical edges of the detector volume are cut away by the requirement 
that the cosmic muon pathlengths be longer than 1 mm. The intersection points are calculated 
with FIND_LINE_REGIDN_INTERSECTION subroutine. 

FIG. 3. Three plastic veto scintillator volumes intersected by tracked cosmic muons. The top 
and bottom stave surfaces, separated by 0.318 cm, are outlined by the cosmic muon intersections. 
The vertical edges of the staves are missing due to the requirement that the cosmic muon path- 
lengths be longer than 1 mm. 

FIG. 4. Main tkoptics X- Window-based input menu prompts a user for a series of detector 
parameters. Predefined geometries are conical, pentagonal, hexagonal and general polygonal de- 
tector shapes. Detector boundary surfaces and wrapping materials could be specified inside this 
window. 

FIG. 5. Block diagram of the optics code logic. The starting photon positions and directions 
are read from an ASCII file or simulated according to a user prescription. The numbers of absorbed, 
"over-bounced" and detected photons are counted. 

FIG. 6. optics simulation of the photon transport in the PIBETA HEX-D1 trapezoidal detec- 
tor. The number of detected photon reflections from the detector surfaces and well as the arrival 
times of the detected photons are shown. 

FIG. 7. optics simulation of the scintillator pulse shape for the PIBETA HEX-D1 trapezoidal 
detector. Only the photon time-of-flight contribution to the pulse-height shape is shown. 

FIG. 8. optics simulation of a two inch photocathode surface illuminated by scintillation 
light. The scintillation photons were generated uniformly throughout the volume of the HEX-D1 
trapezoidal detector, see Fig. [I] and Table |. 

FIG. 9. The simulated (i) and measured (ii) light collection probabilities are shown as functions 
of axial (z) and transverse (x) coordinates for an ideal hexagonal detector (HEX-A) with a two-layer 
Teflon wrapping. The size of the histogrammed two-dimensional bins is 1 x 1 cm 2 . The average light 
collection probability with a 78mmc/> inch phototube is 23%. 
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FIG. 10. (i) Calculated Monte Carlo light collection probability as a two dimensional function 
of axial and transverse coordinates for an ideal half-hexagonal detector (PIBETA HEX-D1) with 
a Teflon wrapping. The average light collection probability with a 46mm</> photocathode is 11%. 
(ii) Second panel shows the measured response of the detector. 



FIG. 11. The two-dimensional light response of the RADPHI lead glass calorimeter monitoring 
plate. The nitrogen laser light is fed into the plate via six 1 mm thick optical fibers at the top and 
bottom sides. 



FIG. 12. (i) Monte Carlo pulse-height spectra of 70MeV photons in the pure Csl modular 
calorimeter. The dashed line histogram represents the case of an ideal detector with uniform light 
collection probabilities. The full lineshape shows the energy spectrum calculated using GEANT code 
and the light collection probabilities calculated with optics program (see Figs. | and [To]), (ii) The 
identical calculation done for 70MeV monoenergetic positrons. 
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optics Simulation I Scintillation Photon Transport 



*a> 
x 

(D 
_C0 

Z5 

CL 

0) 
CD 
O 

> 
< 



2 




2 
4 
6 
-8 
-10 
-12 



-2 



2 



4 



ID 


3001 


Entries 


110 


Mean 


2.960 


RMS 


2.937 


Ao{-exp( 


-t(ns)/3.0) 


10 12 


14 16 



Time (ns) 



FIGURE 7 



32 
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S021-1a: 2D Light Output Uniformity 
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S202- 1 a : 2D Light Output Uniformity 
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